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"The next question was - what makes planets go 
around the sun? At the time of Kepler some peo- 
ple answered this problem by saying that there 
were angels behind them beating their wings and 
pushing the planets around an orbit. As you wih 
see, the answer is not very far from the truth. The 
only difference is that the angels sit in a different 
direction and their wings push inward." 



"Nast§pnym pytaniem bylo - co sprawia, ze plan- 
ety kr§iz§ wokol slohca? W czasach Keplera 
niektorzy ludzie odpowiadali, ze za planetami s§, 
aniolowie, ktorzy popychajg, je machaj§,c skrzy- 
dlami. Jak widzisz, odpowiedz ta nie jest daleka 
od prawdy. Jedyna roznica jest taka, ze aniolowie 
zasiadaj§, w innej pozycji i ich skrzydla pchajg, do 
srodka." 



Richard Feynman (1918 - 1988), 
Character of Physical Law 



Abstract 



The work was carried out in the framework of the KLOE collaboration studying the decays of 
the (f) meson produced in the DA$NE accelerator in the collisions of electron and positron. 

The main aim of this thesis was investigation of the influence of the merging and splitting of 
clusters in decays with the high multiplicity of 7 quanta, which are at most biased by these effects. 
For this aim we implemented the full geometry and realistic material composition of the barrel 
electromagnetic calorimeter in FLUKA package. The prepared Monte Carlo based simulation pro- 
gram permits to achieve a fast generation of the detector response separately for each interested 
reaction. The program was used to study the reconstruction efliciency with the KLOE clustering 
algorithm as a function of the photocathode quantum efficiency. 

It was also used to investigate merging and splitting probabilities as a function of the quan- 
tum efficiency. The conducted studies indicated that the increase of quantum efficiency does not 
improve significantly the identification of clusters. The infiuence of these effects was estimated for 
77 meson decays into Stt" and Kshort meson into 27r". 



Streszczenie 



Praca zostala wykonana w ramach kolaboracji KLOE, ktora zajmuje si§ badaniem rozpadow 
mezonu 0, produkowanego w kolizjach elektronu i pozytonu na akceleratorze DA$NE. 

Glownym celem pracy bylo oszacowanie wplywu efektow "merging" i "splitting" na mozliwosci 
rekonstrukcji rozpadow mezonu 4> w kanaly z fotonami. Wlasnie rekonstrukcja topologii reakcji 
z duz§, zawartosci§t fotonow w kanale wyjsciowym jest najbardziej obci§,zona przez te efekty. W 
tym celu zaimplementowalismy cal§t geometri§ kalorymetru elektromagnetycznego detektora KLOE 
wraz z realistycznym opisem materialow do programu symulacyjnego opartego na metodzie Monte 
Carlo. Program symulacyjny pozwala na studiowanie wybranych kanalow rozpadu oddzielnie. 
Rowniez wysymulowalismy odpowiedz fotopowielaczy na calym kalorymetrze detektora KLOE. 

Zbadane zostaly takze efekty "merging" i "splitting" w zaleznosci od wydajnosci kwantowej 
w fotopowielaczach. Powyzsze studia pokazaly, ze zwi§kszenie wydajnosci kwantowej fotopowielaczy 
nie wplywa znaczg^co na redukcj§ tych efektow. Wplyw efektow "merging" i "splitting" zostal osza- 
cowany dla reakcji rozpadu mezonu 77 na trzy neutralne piony oraz kaonu T^sUort na dwa neutralne 
piony. 
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1. Introduction 



Nowadays, elementary particles and interactions among them are described in the framework 
of the Quantum Mechanics and Standard Model. Tests of these theories and of fundamental sym- 
metries Hke C, P, T and their combinations are therefore crucially important for understanding the 
phenomena in the world of particles [1] . Experiments aiming at the determination of the particle 
properties and search for processes beyond the applicability of the domain of the Standard Model 
are conducted in many particle physics laboratories. 

Such studies are carried out also in the Laboratori NazionaH Di Frascati (LNF) by means of 
the electron-positron collider DA<i>NE and the KLOE detector setup. The KLOE group has taken 
data at the colliding electron-positron center-of-mass energy corresponding to the 4> meson mass. 
The main objectives of the KLOE experimental program are: investigations of the decays of kaons 
produced in pairs from the decays of the <^ meson, the radiative decays of the (f) meson, the rare 
branching ratios of the r] meson and the light mesons particle properties. One of the interesting 
example for the pseudoscalar mesons decay is the G-symmetry violating 77 7r°77 process whose 
branching ratio cannot be described neither in the framework of the Chiral Perturbation Theory 
nor by the Vector Meson Dominance Model. 

At present KLOE detector is being upgrated in order to improve the possibility of the recon- 
struction of the decays of the Ks mesons. In particular in order to investigate a CP-violating 
Ks SttO decay [2]. 

The 77 7r°77 and Ks 37r° are examples of highly interested channels which are however 
very challenging to study experimentally due to the large number of the 7 quanta in the final state. 
This fact is caused by difficulties of the rejection of background channels due to merging and split- 
ting of clusters [3]. For example, for the first mentioned reaction a background channel is 77 ^ Stt'' 
with two merged clusters leading to the same topology as in the case of the rj — ^ 71^77 47. And 
in the case of the Kg Btt^ the background constitues a Ks 27r° reaction with two splitted 
clusters. 

In this diploma thesis we investigate the infiuence of the merging and splitting of clusters on de- 
cays with the high multiplicity of 7 quanta, which are at most biased by these effects. For this aim 
we implemented the full geometry and reaHstic material composition of the barrel electromagnetic 
calorimeter into FLUKA package. The prepared Monte Carlo based simulation program permits 
to achieve a fast generation of the detector response separately for each interested reaction. 

In chapter 2 we will describe the DA$NE collider and the KLOE detector. 

Chapter 3 comprises description of the vertex generator, used to simulate the kinematics of 
physical decays, and the description of the physics models used by the FLUKA Monte-Carlo pro- 
gram for the generation of nuclear, hadronic and electromagnetic reactions. This chapter includes 
also description of implementation of materials composition and geometry of the barrel calorimeter 
in the FLUKA Monte Carlo. 

Chapter 4 presents a calibration of the DIGICLU program used for the reconstruction of the 
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photomultiplier response and for the cluster recognition. This chapter includes: i) preparation 
of the data sample for e+e^ ^ 7^7 ^ 87 reaction, ii) estimation of the attenuation length 
for scintillating fibers, iii) caHbration of ionization deposits and implementation of the threshold 
formula to the source code. 

In chapter 5 we describe effects of merging of clusters for 77 Stt" as a background for the 
rj — *■ 7r°77 channel and splitting of clusters for Ks — > 27r° as a background for Ks 37r° reac- 
tion. Further on reconstruction efficiency with the KLOE clustering algorithm as a function of 
the photocathode quantum efficiency is presented. Finally, merging and splitting probabilities as 
a function of the quantum efficiency are studied. 

Chapter 6 summarises the whole thesis and brings the conclusions and remarks. 

This thesis is supplemented with appendics where section A presents a kinematic fit procedure 
and section B generally describes Monte Carlo Methods. In appendix C energy distribution for 7 
is presented. Section D presents a time distribution for single and multi-gamma hits. The sections 
E and F contain an estimation of probability for multi-gamma hits at a single calorimeter module 
and example of the event reconstruction for the process e+e~ ^ ^ 777 ^ 37r°7 — *■ 77, respec- 
tively. The G section comprises a description of the energy deposition in the barrel calorimeter as a 
function of azimuthal angle and the last section presents definition of the used coordinate system. 



2. KLOE2 at DA$NE experimental facility 



The KLOE (Klong Experiment) detector is installed at the interaction point of the electron 
and positron beams of the DA$NE (Double Annular (/)-factory for Nice Experiments) collider 
operating in the Laboratori Nazionali di Frascati (LNF). 

It has been designed with the primary goal to measure the CP violation parameter R(y) [4] 
with a sensitivity of one part in ten thousand by using the double ratio method [5] . 

This detector was fully constructed by the end of the year 1998 [1], and since then it was taking 
data for seven years. The experimental program was completed with integrated luminosity of 
2.5 fb~^ obtained predominantly with the center-of-mass energy equal to the mass of the (j) meson 
(Vs M4, = 1019.456±0.020 MeV [6]). 



Decay 


BR(%) 


(j) K+K- 


49.1 


4> ^ KsKl 


33.8 




15.6 


?77 


1.26 



Table 2.1: Main decays of the (p meson [59]. 

The cross section for the production of the 4> vector meson is large and amounts to 
cr(e+e~ ^ 4>) = 3.1 ph. 

The (j) meson decays predominantly into pairs of neutral or charged kaons (see Table 12. ip . 

2.1 DA$NE collider 

DA$NE at Frascati consists of two intersecting crossing accelerator rings, one for positrons 
and one for electrons. The main DA$NE parameters are presented in Table [221 and scheme of 
the facility is shown in Fig. 12.11 



Parameter 


Value 


Energy [GeV] 


0.51 


Trajectory length [m] 


97.69 


RF frequency [MHz] 


368.26 


Bunch length [cm] 


1-3 


Number of colliding bunches 


111 


Emittance, Cx [mm-mrad] 


0.34 



Table 2.2: The main characteristics of the DA$NE collider ]7]. 
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Positrons and electrons are accelerated in the LINAC (Linear Accelerator) which delivers 
electron or positron beams in the energy range from 25 to 725 MeV with intensities varying from 
10^" particle per pulse down to a single-electron [7]. This 60 meters long accelerator is the heart 
of the DA$NE injection system. It is an S-band machine (2.865 GHz) which delivers 10 ns pulses 
at a repetition rate of 50 Hz. Electrons, after acceleration to final energy in the LINAC, are 
accumulated and cooled in the accumulator and transferred to a single bunch into ring. Positrons 
require first accelerating of electrons to about 250 MeV to target in the LINAC, where positrons are 
created. Afterwards the positrons follow the same accelerator elements as electrons [1]. Positrons 
and electrons after acceleration and accumulation process run around in two storage rings and hit 
in the collision points. This facility is called a Frascati ^-Factory complex because it was able to 
produce about 3.3 bilious of mesons during the years from 2000 to 2005. The KLOE experiment 
was located in one of the two collision points at DA<I>NE collider, whereas the second collision 
point was alternatively occupied by two other experiments: DEAR [8] and FINUDA [9]. 

At present a new e+e" interaction region is being constructed [10] in order to increase the 
collider luminosity up to an order of magnitude, from 10'^^cm~^s~^ to « 10'^^cm~^s~^ [3]. 




Figure 2.1: Scheme of the DA$NE collider. The figure is adapted from [1]. 

At the interaction point (IP) the beam pipe has the shape of a sphere which is made of a 
beryllium-aluminium alloy with 10 cm diameter and 50 thickness. The beryllium, having a low 
atomic number, has been used to minimize the interaction of particles produced at the interaction 
point with the beam pipe material [6]. 

2.2 KLOE detector 

The KLOE detector, which is shown schematically in Fig. 12. 2^ was designed for the study of 
CP violation in the neutral-kaon system. It consists of two main elements: i) an electromagnetic 
calorimeter (EmC) [11] for the detection of 7 quanta, charged pions and Kl mesons and ii) a 
large drift chamber (DC) [12] for the measurement of charged particles trajectories [13]. The drift 
chamber and the calorimeter are inserted in the field of the superconducting coil which produces 
a magnetic field parallel to the beam axis [11]. The field intensity is equal to 0.52 T [6]. 
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Figure 2.2: The KLOE detector. For the description see text. The figure is adapted from [3]. 

The radius of the active part of the KLOE detector is two meters. This size enables to register 
about 40% decays of neutral long-lived kaons [1]. In the near future KLOE will be upgraded by 
inner tracker and 77 tagger in order to improve its tracking capabilities [3]. 

2.2.1 Drift Chamber 

The KLOE Drift Chamber (DC) [12] consists of 12582 drift cells (2x2, 3x3 cm^) arranged in 58 
cylindrical layers surrounding the beam pipe. The diameter and length of the DC is equal to 4 m 
and 3.3 m, respectively [14]. It is filled with 90% hehum and 10% isobutane gas mixture, giving 
a radiation length (gas + wires) equal to 900 m. Charged particles travehng through the drift 
chamber are ionizing gas medium and then electrons created along the particle trajectory drift to 
the wires with positive voltage. A multiplication mechanism causes detectable signal at the wire's 
end [1]. The DC measures the Kl,s charged vertex with ^ 1 mm accuracy. It provides fractional 
momentum resolution of ~ 0.5% for low momentum tracks. It is transparent to 7 down to 20 
MeV and limits to acceptable levels the regeneration and multiple scattering. 




Figure 2.3: Photo of the KLOE drift chamber. The figure is adapted from [15]. 
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The hit position resolution is 150 in the central part of the cell, increasing close to the wire 
and towards the cell boundary. The spatial resolution in the "0 coordinate" (azimuthal angle) is 
well below 200 ^m. The hit identification efficiency is larger than 99%, whereas the efficiency for 
associating an existing hit to track amounts to about 97% [14]. 

2.2.2 Electromagnetic Calorimeter 

The KLOE calorimeter is made of the lead layers with about 1,2 mm thickness (200 layers per 
1 module) which are filled with scintillating fibers of 1 mm diameter. The whole electromagnetic 
calorimeter (EmC) consists of three main parts: barrel and two endcaps. Barrel (Fig. 12. 4p is 
composed of 24 modules with trapezoid shape of 23 cm thickness aligned with the beams and 
surrounding the drift chamber detector. Endcaps are situated over the magnet pole pieces (see 
Fig. 12. 2p and hermetically close the calorimeter with 98% of 47r [1]. Each of the two endcaps 
calorimeters consists of 32 vertical modules with length ranging from 0.7 to 3.9 meters. The 
endcap modules are bent and their cross-section with plane parallel to the beam axis is rectangular 
with thickness of 23 cm [16]. 




Figure 2.4: Photo of the KLOE calorimeter. One sees 24 modules of the barrel and the inner plane of 
the endcap. The figure is adapted from [15]. 

The volume of the calorimeter consists of 50% fiber, 40% lead and 10% of glue. The measured 
performances for this detector are: full efficiency for 7 quanta from 20 to 500 MeV [14], a{x) ^ 1 cm, 

2.2.3 Upgrade of the KLOE detector 

At present the KLOE detector is being upgrated in view of the new experimental program which 
will extend the studies to the more precise measurement of the Ks mesons and the production 
of meson in the 77 fusion [3]. For this aim a vertex detector and the 77-tagger [17] are being 
built [26]. 

Vertex Detector 

In the KLOE detector the first hit was measured by drift chamber at a radius of 28 cm from 
the interaction point (IP) [3] . Therefore in order to improve the resolution of the determination of 
the Ks and Kl decay near the interaction point a new vertex detector is constructed [17]. This 
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detector consists of five concentric layers of cylindrical triple-GEM (C-GEM), completely realized 
with very thin polyimide foils [18-22]. The scheme of the vertex detector is presented in Fig. 12.51 

70 cm 




Figure 2.5: Schematic view of the vertex detector. The figure is adapted from [3]. 

It will be made of light materials in order to minimize 7 absorbtion. The resolution on r-4> 
surface is expected to be around 200 ^m and on z coordinate about 500 /im [23]. 

77 tagger 

The main task of this detector is the detection of e+ and e~ from 77 reactions emitted at 
small angle [23,43] with the widest possible energy ranges. 77 tagger will provide information on 
the angle and the energy of the scattered electrons and positrons [24] and hence it will permit to 
study the production of mesons in 77 fusion via the reaction: 

e+e" ^ e+e^7*7* e+e" + X, (2.1) 

where X is some arbitrary final state allowed by conservations laws. The tagger consists of mi- 
crostrip silicon detector and plastic scintillator hodoscope. 



New QCAL 

The upgrade of QCAL detector (Fig. 12. 6|) [25] was needed, because the interaction region was 
modified and in the present scheme the angle between colHding beams has been increased from 8 
to 18 degrees, which practically excludes the possibility to use the existing QCAL calorimeter [26]. 
The schematic view of upper part of the previous QCAL detector is presented in Fig. 12.61 
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Pb layers {1 .9 mm) 




Figure 2.6: Schematic view of the QCAL detector. The figure is adapted from [25]. 

The upgraded detector, with improved position resolution, will allow to extend the search for 
Kl 27r° events also in case when three photons are reconstructed in the EmC and one photon 
in QCALT, thus strongly reducing the correction for acceptance [26] and the contamination from 
Kl Stt" when some photons hit the detector. 



3. Simulation tools 



In this chapter we will describe tools used for simulations of signals in the EmC. To compute 
the physical response of the barrel calorimeter we applied the following programs: i) VERTEX 
GENERATOR which simulates kinematics of the physical reactions ii) FLUKA package which 
reproduces particle interaction, propagation and a realistic light output in the scintillating fibers, 
and iii) DIGICLU reconstruction program which simulates the response of photomultipliers and 
reconstructs energy and time for particles hitting the calorimeter. FLUKA package is mostly used 
for low energy physic, like: hadronic interactions, medical science, neutrino beam simulation or in 
dosimetry science. Using FLUKA we can achieve accurate description of the physics processes in 
the KLOE detector. 

3.1 Vertex generator 

In order to enable a fast simulations of the meson decays in the KLOE detector we have 
prepared a "vertex generator" for simulations of the final state momenta of the decay productl^. 
This program is compatible with FLUKA Monte Carlo and simulates decays of 4> meson to channels 
with ry, 77' and kaons. It is based on the CERN library procedure GENBOD which generates four- 
momentum vectors of particles in the final state which are homogenously distributed in the phase 
space [27]. This tool permits to calculate four- momenta of final state particles in the rest frame of 
the decaying object. Additionally we implemented an angular distribution for 7 quanta 

from the radiative <j) meson decay (0 is a vector meson with spin = 1). At present, for the purpose 
of this work the generator can simulate six below Hsted reactions: 



1. 


e+e~ ( 




37, 




2. 


e^e~ ( 




37r07 - 


.77 


3. 


e^e^ —t ( 




.37 , 




4. 


e^e^ < 




> 37r07 - 


-.77 


5. 


e+e" ( 


i> ^ KsK 


L , 




6. 




f> 7r'^7 - 


-.37 . 





It is however constructed in a way which enables an easy implementation of any further process. 

3.2 FLUKA - a multi-particle Monte Carlo transport code 

This section is based on the informations given in references [29,32]. FLUKA is a multipurpose 
transport Monte Carlo code, able to treat hadron-hadron, hadron-nucleus, neutrino, electromag- 
netic, and ^ interactions up to 10000 TeV [28]. Also the transport of the charged particles is well 

^In our group this generator is called jjpluto. 
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recreated (handled in magnetic field too). This Monte Carlo is based, as far as possible, on orginal 
and well tested microscopic models. Due to this "microscopic" approach to hadronic interaction 
modeUing, each step is self-consistent and has well confirmed physical bases [28]. FLUKA can be 
used in realistic simulations involving following physical effects: 

Hadron inelastic nuclear interactions - are based on resonance production and decay be- 
low a few GeV, and above on the Dual Parton model. FLUKA can also simulate photo-nuclear 
interactions (described by Vector Meson Dominance, Delta Resonance, Quasi-Deuteron and Giant 
Dipole Resonance models) [29]. 

Elastic Scattering - this process is accomplished with tabulated nucleon-nucleus cross sec- 
tions. Tabulated are also phase shift data for pion-proton and for kaon-proton scattering [29]. To 
this Monte Carlo also parametrised nucleon-nucleon cross-sections were implemented. 

Nucleus-Nucleus interactions - are treated through interfaces to external event generators. 
From 0.1 to 5 GeV per nucleon a modified RQMD is used, whereas for higher energies DPMJET 
generator is applied [30]. 

Transport of charged hadrons and muons - a code can handle electron backscattering and 
energy deposition in thin layers even in the few keV energy range. The energy loss is simulated 
with Bethe-Bloch theory with account for spin effects and ionisation fiuctuations [29]. 

Low-energy neutrons - for neutrons with energy lower than 20 MeV, FLUKA uses neutron 
cross-section Hbrary, containing more than 140 different materials. The transport of these parti- 
cles is realised by simulation of the standard multigroup transport with photon and fission neutron 
generation [29]. Detailed kinematics of elastic scattering on hydrogen nuclei is implemented too. 

Electrons - FLUKA uses an original transport algorithm for charged particles, including a 
complete multiple Coulomb scattering treatment giving the correct lateral displacement even near 
boundary. The Landau-Pomeranchuk-Migdal suppression effect and the Ter-Mikaelyan polarisa- 
tion effect in the soft part of the bremsstrahlung spectrum are also implemented. Electrons are 
propagated taking into account a positron annihilation in fiight and Moller scattering effects. 

Photons - to reproduce photons physics the following effects were implemented to this Monte 
Carlo: pair production with actual angular distribution of electrons and positrons, Rayleight 
scattering, photon polarisation effect, photo-hadron production, Compton effect with account for 
atomic bonds through use of inelastic Hartree-Fock form factors. 

Optical photons - generation and transport is based on Cherenkov, Scintillation and Transi- 
tion Radiation. 

Neutrinos and muons - FLUKA package includes also cosmic ray physics [31] and can sim- 
ulate muons and neutrinos interactions. This tool is used for basic research and applied studies in 
space and atmospheric fiight dosimetry and radiation damage [31]. Interactions are implemented 
independly of tracking procedures [29] . By default they are tracked without interactions [29] . 

3.3 Implementation of the calorimeter material composition in 
FLUKA 

The Monte Carlo code FLUKA is used to determine the position, time and energy of the 
ionization deposits in the fibers caused by particles hitting the calorimeter. We used the FLUKA 
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"lattice" tool to design the fiber structure of the calorimeter module [32]. In the base module the 
calorimeter is simulated in detail (see Fig. IS.ip , both under the geometrical point of view and with 
respect to the used materials. 

LEAD 




replicas 




200 layers 

Figure 3.1: Details of the implementation of the calorimeter layers [33]. 

The calorimeter is built out of lead, plexi and glue [34] which are composed as follows: 

• lead compound is made of Pb (95%) and Bi (5%), 

• the scintillating fibers are made of Polystyrene {C2H^), 

• the glue is a mixture of the epoxy resin {C4H4O, p = 1.14 g/cm^) and hardener 
(p = 0.95 g/cm^) constituting 72% and 28% of the mixture, respectively. 

The hardener has been simulated as a mixture of several materials which composition is given in 
Table [331 



Compound 


formula 


fraction 


Polyoxypropylediamine 


C7H20O3 


90% 


Triethanolamine 




7% 


Aminoethylpiperazine 




1.5% 


Diethylenediamine 


C4H10N2 


1.5% 



Table 3.1: Composition of the hardener mixture [34]. 



The module consists of two types of scintillating fibers: "Kuraray" and "Pol. hi. tech". The 
first are implemented until a depth of 12 cm and the remaining part of the module is built with 
"Pol. hi. tech" fibers [34]. A proper attenuation lenght parameter for fibers material in order to 
calculate a light intensity at each side of the module was taken into account. The attenuation 
function of fibers was described with the following formula [34] : 

B ^ A- + {1 - A) ■ , (3.1) 

where y is distance between the place of deposited energy and the photocathode. The attenuation 
factor B is the ratio between detected and generated light signal which changes as a function of the 
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distance (on y axis) from the generation point to the photocatode. The value of the parameters 
for "Kuraray" and "Pol. hi. tech" are shown in Table [331 

However, it is worth mentioning that, as we will show in the next chapters the data can be also 
well described using this formula but with only one exponential function. 





A 


Ai cm 


A2 cm 


Kuraray 


0.35 


50 


430 


Pol. hi. tech 


0.35 


50 


330 



Table 3.2: Fibers parameters. 



The energy deposits are computed by FLUKA [35] taking into account the Birks effect [36] 
(see equation 13. 2p , that is the saturation of the Hght output (L) of a scintillating material when 
the energy release is high [35]. For high densities of energy deposition, due to the quenching 
interactions between the excited molecules along the path of incident particles, the light output is 
not changing linearly with the energy deposition but instead it can be described as [37-40]: 



dL 

dx 



ax 



(3.2) 



This law describes the Hght output of (organic) scintillators [39,40]. The "k" constant depends on 
the particle type, and is in the order of 0.01 g cm^ MeV~^, and the parameters ci and C2 are equal 
to [33,41,42]: 



ci = 0.013 cm^ MeV"^ , 
C2 = 9.6 X 10"^ cm^ MeV"^ 



3.4 Implementation of the whole geometry of the barrel calorimeter in 
FLUKA 

FLUKA simulation package is mostly used in low energy physics, for example in: medical 
science, dose evaluation, neutriono beam simulation and for simulations of the low energy hadronic 
interactions. Therefore, the program is especially suited for simulations of the response of the 
KLOE electromagnetic calorimeter, because at KLOE experiment we deal with particles in the 
rather low energy range up to 1 GeV. 

The main idea for implementation of the geometry of the KLOE detector with FLUKA was 
proposed by Giuseppe Battistoni [44]. It was first realized for one calorimeter module with 
rectangular shapjfl, with accurate description of materials. The base cell with dimensions of 
52 cm X 1.2 mm x 430 cm consists of the lead block filled with 385 scintillating fibers and glue 
cylinders. This cell was replicated 199 times (using FLUKA lattice tool) to build 200 layers forming 
the calorimeter module (see left panel in Fig. 13. 2p . 

Our first step in the way to implement the whole barrel of electromagnetic calorimeter into 

^Implementation to FLUKA package was done by: G. Battistoni, B. Di Micco, A. Ferrari, A. Passeri and 
V. Patera. 
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FLUKA was the extension of the rectangular module to the realistic trapezoid shape. This was 
achieved by declaration of a new base cell consisting of a scintillating fiber and glue cylinder both 
inserted into a small block of lead. The length of the cell is equal to the length of the module 
(4.3 m) and its cross section corresponds to the square of 1.2 mm, where the fiber diameter is equal 
to 1 mm and width of the glue cylinder equals to 0.1 mm. Then this new base cell was replicated 
about 4000 times to build two triangular sections on the left and right side of the module. A 
complete visualisation of the new trapezoid geometrjfl is shown in Fig. 13.21 

RECTANGULAR GEOMETRY TRAPEZOID GEOMETRY 





E 



X rcml X Fcnil 

Figure 3.2: Cross section of the EmC module with rectangular and trapezoid geometry. Visualization by 
the FLAIR program [45] using as an input a geometry setup files of FLUKA. 

In order to raise our confidence to the functioning of the programs we have simulated energy 
deposits in a single calorimeter module for the reaction: e+e~ (f) ^ ^ Sn^-f 7-f. The 
results of these simulations for 20000 events are shown in Fig. | 
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Figure 3.3: Distribution of energy deposits in scintillating fibers of a single unit of the electromagnetic 
calorimeter. 

^Visualization is based on the FLUKA sources and prepared using a FLAIR (FLUKA advanced visualization 
geometry tool) [45]. 
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A trapezoid structure which was built with a realistic light output from the scintillating fibers 
can be easily recognized. One can also see that as expected the most energy deposits are observed 
in the bottom part of the module. As a next step, it would be natural to replicate 23 times the 
defined geometry for one module and to build a full barrel of the electromagnetic calorimeter. This 
task was unfortunatelly impossible to realize with present version of FLUKiJ^ [32] because this 
version doesn't permit to replicate the region which had been replicated before. It is due to the 
fact that the lattice repHcation on the second and higher levels is not implemented yet. This is the 
reason why geometry was built in another, unfortunatelly much more complicated way. 

The barrel was first defined as 24 empty volumes (containers). In one of these volumes we 
implemented the base cell which described one layer of lead, fibers and glue. This cell was replicated 
199 times in this volume, to build a rectangular section part of the main trapezoid module (on 
the top in Fig. I3.4|PI . and these base cell was repHcated 200 times in each of the remaining 23 
modules. Next in the first module we filled the full structure of lead, fibers and glue cylinder in 
the two triangle areas. These two regions were then replicated at the corresponding positions in 
the remaining modules [46]. 

Fig. 13.41 shows a visualisation of the whole barrel calorimeter geometry made with FLAIR [45]. 
The left panel of the next figure presents details of the edge area, between base module on the top 
and the first replicated area on the right side of it. In the right panel of Fig. l3.5l the energy deposits 
in scintillating fibers for this area are presented for e+e^ cj) ^ rj-f ^ 3n^-f 7j reaction with 
statistics of 1000 events. 
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Figure 3.4: Visualization of the 24 trapezoid modules of the barrel calorimeter with FLAIR. 

Using this geometry we are able to study in details energy deposits also at the edges of the 
modules of the calorimeter, where there are small areas without fibers [47,48]. 
"FLUKA 2006 for GNU linux operating system. 

^The same idea was used to build a single module with rectangular shape. 
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Figure 3.6: Energy depositions in fibers in the KLOE barrel calorimeter. 

Finally in Fig. 13.61 as an example of the proper implementation of the geometry we show energy 
deposits in scintillating fibers in the whole barrel calorimeter. 

3.5 Simulations of the photomultipliers response 

After simulation of energy deposits in scintillating fibers with FLUKA program we calculated 
amplitude of signals in photomultipliers using DIGICLU program. This program also enables 
to reconstruct the time, position and energy for all readout segments (cells) of each calorimeter 
module [11]. Each module is divided into 60 cells which are situated in 5 layers and 12 columns 
(see Fig. 13. 7p . Each cell is read out on both sides by photomultiplers. This segmentation provides 
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the determination of the position of energy deposits in the calorimeter. Altogether, the barrel 
calorimeter consists of 1440 cells, which are read out by photomultiplayers on two sides (referred 
to as side A and side B in the following). 
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Figure 3.7: Schematic view of the readout cells structure on the one side of the barrel module. Filled 
circles represent photomultipliers. 

For each cell two time signals T"^, (digitized by the Time to Digital Converter (TDC)) and 
two amplitude signals 5*^, 5"^ (measured by Analog to Digital Converter (ADC)) are recorded. 
The arrival time t and position s of the impact point along the fiber direction (the zero being taken 
at the fiber center) is calculated with the aid of times measured at two ends as: 

tins) = l{t^ + t^ -t^-t^)-L. , (3.3) 
sicm) = ^(t^-t^-t^ + ti) , (3.4) 

with t^'^ = c"''' • T^'^ , where c"^'^ are the TDC calibration constants, t^'^ denotes overall time 
offsets, L stands for length of the cell (cm) and v is the light velocity in fibers (cm/ns). The energy 
signal is calculated according to the formula: 

gA,B _ gA,B 

Ef'^ (MeV) = ^— ^ • K, (3.5) 

where 5*0, i are the zero-offsets of the amplitude scale, Sma corresponds to the response for the 
minimum ionizing particle crossing the calorimeter center and K factor gives the energy scale in 
MeV, and it is obtained from signals of particles with known energy. 

The total energy deposited in a cell is calculated as the mean of values determined at both 
ends for each cell [11]. The energy read-out has been simulated by including both the genera- 
tion of photoelectrons with the Poisson distribution and the threshold on the constant fraction 
discriminator. 

3.6 Description of the KLOE clustering algorithm 

As we presented in the previous section the DIGICLU program simulates photomultipHers 
response and reconstructs energy deposit, position and time for particles passing through each cell. 
These values are used to recognize groups of cells (clusters) belonging to particles entering the 
calorimeter. For this aim a clustering algorithm is used [11]. Ideally, to each particle it should 
assign exactly one cluster but in practice it is not always the case. After the recognision of clusters 
the program reconstructs the spatial coordinates and time of each shower with high accuracy, 
needed to reconstruct the decay vertex of the Kl [49]. In particular the algorithm is based on the 
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following steps. First for each cell the position and energy of the shower is reconstructed. Next 
preclusters are built by connecting the neighboring cells in time and space in order to recreate a 
full shower [11]. Cells are taken into account in searching of preclusters only if times and energy 
signals are available on both sides, otherwise these cells in the most cases are added to the already 
recognized clusters [6]. Subsequently, preclusters are splitted if the spread of the time of the 
assigned cells is larger than 2.5 ns. On the other hand cells are merged in one cluster if a distance 
between them and the center of the precluster is less than 20 cm. After this check the groups of 
cells are defined as clusters which position and time are computed as energy-weighted averages of 
the contributing cells. 
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4. Adjustment of the detector properties in simulations 

In order to simulate a realistic response of photomultipliers we have adjusted parameters used 
in the programme by comparing the output of simulations to results from the experimental sample 
of events identified as <j> ^ f]j ^ 3j process. 

4.1 Preparation of the data sample {e~^e~ — > ^ 777 ^ 87) 

At a first stage, an experimental data sample for the e+e^ ^ 87 reaction has been 
extracted [50] applying the kinematic fit procedure to the 3 7 events with the following conditions: 

3 

3 
i=l 

U^^ = 0, 

(4.1) 

where, c is the light velocity, ti denotes the time of each reconstructed cluster and r^ stands for 
the distance from the vertex collision point to the cluster centroid position. Ecj, and denote the 
total energy and momentum of the (j) meson, respectively. The first and second condition results 
from energy and momentum conservation rules. The third requirement ensures that we take events 
which orginate only from the collision point. 




Figure 4.1: The distribution of the square of the invariant mass for 7 pairs from e^e ^ — > 87 reaction 
from the experiment (left panel [50]) and from the simulation (right panel). 
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As a next step, the distributions of the square of the invariant mass for pairs of the 7 quanta 
for e+e^ ^ ^ 87 process have been constructed. The results are shown in the left panel of 
Fig. 14.11 The numbers were assigned to the 7 quanta such that they are ordered according to the 
increasing energy: Ei < E2 < E3. Thus 73 is the photon with the highest energy and 71 with the 
lowest. The right panel presents results of simulations for the following reactions: 

1. e+e^ ^ ^ 7/7 ^ 87 , 

2. e"'"e" (j) ^ 7r°7 — > 87 , 

3. e+e" — ^ 87 . 

The blue bevel line and the vertical band at the ?/ mass squared denote signals from the first 
reaction. The channel where <j) meson decays to 7r°7 is described by the red line with value of 
m^7i72 equal to 0.018 GeV^. And finally the green "triangular" region orginates from reactions 
where 4> meson decays directly to 87 quanta. 

For the further analysis we took all events from the vertical 77 band which are positioned below 
the red bevel fine (see left panel in Fig. 14. ip [50]. 

4.2 Determination of the simulation parameters 

In this section we describe determination of the attenuation length for the scintillating fibers, 
determination of the threshold function, and caHbration of the ionization deposits in the electro- 
magnetic calorimeter. 

4.2.1 Estimation of the attenuation length of the scintillating fibers 

The first step for estimation of the attenuation length was the selection of events corresponding 
to the — > 777 — > 777 reaction and identification of monoenergetic 7 from the (/) ^ 777 radiative 
decay. In the next step we simulated the total energy deposited in cells of a given plane as a 
function of Z coordinate [51]. 

We used the following selection cuts: 
1. Monochromatic photons (361 < E^ < 365 MeV) from (j) meson decay , 

2 I 2 reconstructed _ ^ generated | ^ 2Q Qm 

The first condition guarantees that we take only events with 7 quanta from (j) meson decay 
directly. The second condition guarantees that we take only events where the distance between 
position where 7 quanta hit the module in simulation and reconstructed position of the cluster 
must be less than 20 cm. 

The results for 5 layers are shown in Fig. 14.21 One can see that distribution of the average energy 
depositions in cells are not uniform, and therefore for the studies of the attenuation length it is 
mandatory to divide the experimental data by the Monte Carlo distributions. 
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Figure 4.2: Simulation of the total energy deposits in cells for each of the 5 calorimeter layers for the 7 
quanta from the ^ 777 decay. Courtasy of E. Czerwinslci and B. Di Micco [54]. 

In case of the first layer, which is the nearest one with respect to the collision point, the shape 
of the energy distribution is mostly related (i) with the angular distributions of the 7 quanta from 
the (j) — *■ 7771 — *■ 717273 decay, (ii) with the changes of the solid angle distribution along the module, 
and (iii) with the change of the input angle into calorimeter surface. The angular distributions for 
these particles are presented in Fig. 14.31 




Figure 4.3: Simulated angular distributions for photons from reaction: e^e — > 777 ^ 87. 
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The angular distribution for 7 quanta from the (p ^ rj^ decay is described with the formula: 
i+eos e ^ Interestingly, Fig. 14.31 indicates that also 7 quanta from subsequent 77 meson decay have 
still small contribution in the form of cos^O. The decrease along the module of the solid angle as 
seen from the interaction point is shown in Fig. 14.41 
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Figure 4.4: The solid angle distribution as a function of the distance from the center of the module. 

The variation of the solid angle decrease the events population towards the edges. The third 
effect which play a role in this context is the angle at which particle hits the detector surface and 
hence the angle of the shower propagation (see Fig. 14. Sp . 



Collision point 



Figure 4.5: Schematic view of the particle tracks for two different entering examples. 

Due to this fact in the first layer the energy deposited by showers in the middle of the module 
is lower than on the sides. Therefore, in subsequent layers relatively larger fraction of the energy 
of the shower will be deposited in the middle of the module and less on its sides. Hence, the shapes 
of distributions observed in Fig. 14.21 is determined by three discussed effects: angular distribution 
of 7 quanta, changes of the solid angle, and variation of the fractional length of the showers in the 
layers as a function of the impact angle. 
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Figure 4.6: The energy distributions as seen at one side of the module: Monte Carlo (left column) and 
experimental data (right column). The curves are described in the text. Courtasy of Eryk Czerwihski and 
B. Di Micco [541. 
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In order to achieve a uniform, energy deposition along the module a simulated energy deposit 
distributions (Fig. 14. 2p were fit with second order polynominal functions which were then used to 
normalize the experimental spectra and for the check also simulated spectra were normalized in 
the same way. Fig. 14.61 presents an energy of signals registered at one side of the module as a 
function of the distance between the edge of the module and the hit position. In order to establish 
a value of the light attenuation length (A) we fit to the spectra from Fig. 14.61 (right column) an 
exponential function: 

Eiy)^a-e-^, (4.2) 

with a and A being a free parameters [54]. The result of the fit is shown as dark blue, continuous 
line in Fig. 14.61 (right column) . The weighted averages of attenuation length determined at both 
sides (A and B) amounts to 522±7 cm and 407±9 cm, for 1,2,3 layers and for 4,5 layers, respectively. 
The result is shown as a red, dashed curve in Fig. 14.61 (right). 

As a next step the determined attenuation length were used to simulate the light output at the 
edges of the module with the DIGICLU program. The result is shown in Fig. 14.61 (left). These 
studies were made to check if results from simulations are in agreement with results which were 
achieved from the experimental data sample. The code of the lines is the same as in the right panel. 
The simulated and experimental spectra are in very good agreement except a region below 50 cm 
and above 400 cm. This is due to the fact that this part of the barrel is situated behind the endcap 
modules. And in the simulations presented in this thesis only a barrel part of the calorimeter 
is taken into account. Furthermore we tested also a fit with a function with two attenuations 
lengths [34]: 

E{z) = ai ■ + a2 ■ e"^ . (4.3) 

The result is presented in Fig. 14.61 where left column indicates results from simulations and right 
column shows experimental data. As one can see both presented fits describe the data well. 

4.2.2 Parametrization of the energy threshold function 

In the experiment the signals which amplitude is lower than a given threshold are not read by 
the TDC board and information about their times is missing. In order to account for this effect in 
simulations the threshold curve was determined using the ratio between the number of cells with 
time information over all cells as a function of the cell energy. This ratio is shown in Fig. 14.71 
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Figure 4.7: Ratio of number of cells with both time and energy signal to total number of cells with an 
energy signal. Courtasy of E. Czerwihski and B. Di Micco [54]. 

One can see that the distribution is not a step-Hke function but rather a smooth increase of 
the ratio in the range from to 2 MeV is observed. The shape of this increase can be described 
by the Fermi-Dirac distribution [52,53,55]: 

ftH{E) = il+e-^r\ (4.4) 

with free parameters fj, and a. A fit for the range of energy near the threshold resulted in the 
following values: 

cr = 0.2644 ±0.0068 MeV , 
/i = 0.5648 ± 0.0087 MeV . 

4.2.3 Calibration of the ionization deposits 

The next step in the tuning of the simulation program is the adjustment of the absolute scale 
of the energetic response of the calorimeter. For this purpose we compared results of simulations 
with experimental distribution of energy deposits in the calorimeter module. The results for five 
layers are shown in Fig. 14.81 (left column). 
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Figure 4.8: Comparision of energy deposits in 5 layers as determined from simulation (black line) and 
experiment (red line). The presented results before (left column) and after (right column) calibration. 
Courtasy of E. Czerwinski and B. Di Micco [541. 



4.2. Determination of the simulation parameters 



39 



One sees that the energy deposits reconstructed for simulated events are much smaller then 
for the experimental data. In order to correct this discrepancy, the simulated values of energy 
deposits were multipHed by the factor of 3.74. The shape of the corrected spectra is in a very good 
agreement with experimental distributions for each layer (Fig. 14.81 right). 

Next, correction was needed because the clustering algorithm not always is able to reconstruct 
the entire energy deposited in the clusters. That's why sum of energy deposits which were simulated 
is greater than sum of energy in reconstructed clusters. For this case a linearity correction was 
needed and it was obtained with requiring straight Hue with slope 1 (blue Hne in Fig. 14. 9p . 
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Figure 4.9: Calibration of the simulated and reconstructed energy. Courtasy of E. Czerwiriski and B. Di 
Micco [54]. 

The red curve represents the situation before caHbration (slope parameter is not equal 1) and 
blue after. Applying the corrections described in this section we achieved very good agreement 
between data and simulations of energy response of the calorimeter when using FLUKA and 
DIGICLU procedures [54]. 
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5. Estimations of background due to the merging and 
splitting of photon clusters 



The main motivation of this work is the estimation of merging and splitting effects for reactions 
Uke e.g. 77 3iT° or Ks 27r° which may constitue background in studies of the 77 tt^'^j and 
Ks 37r° reactions, respectively. In Table lOl we present a Branching Ratios (BR) for the main 
investigated reactions and for the background channels. 



Main Decay 


BR (%) 


Background Channel 


BR (%) 


Effect 


rj —> 7r'^77 


(4.4 ± 1.6) •10"'^ 


r;-> 37r° 


(32.51 ± 0.28) 


merging 


Ks BttO 


< 1.2 -10-9 


Ks 27r" 


(30.69 ± 0.05) 


splitting 


Kl ^ 2^0 


(8.69 ± 0.04) •lO"^ 


Kl Stt" 


(19.56 ± 0.14) 


merging 



Table 5.1: Branching Ratios for examples of studied reactions and the background decays [59]. 



In this Table one can see that all the background reactions have much bigger BR (many orders 
of magnitude) than the investigated channels. Therefore, even small merging and splitting effects 
may obscure considerably observation of signals from 77 7r°77, Ks — + Btt" or Kl '2n° decays. 

5.1 Merging of clusters from r] Stt^ as a background for the rj n^-jj 
channel 

The example of the merging effect, inducing a wrong cluster counting, comes from analysis of 
the ?7 7r"77 events [6]. The Branching Ratio (BR) for 77 3n^ is by several orders of magnitude 
more frequent than 77 — > 7r°77 (see Table \5A} . 
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Figure 5.1: Invariant mass of four photons for the rj tt 77 events and for the 77 — > 3n background 
with two merged clusters [6]. 
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5.1. Merging of clusters from rj Stt" as a background for the rj — > 7r"77 channel 42 

In Fig. 15.11 a consequence of the merging of clusters is presented for the case when 67 quanta 
from 37r° decay were reconstructed as 4 clusters. In this case the topology of 77 37r° becomes 
equal to that of the i] — ^ tt^j'j reaction. And as a result on invariant mass distribution one can 
see a big peak from background channel exactly at the mass of the r] meson [3]. 

The next very interesting example of reactions for which problems with merging of clusters 
is significant and which can be studied with the KLOE calorimeter [11] is CP violating decay 
Kl Stt" . An example of such event is presented in Fig. 15.21 




Figure 5.2: A CP-violating KsKl event. The Ks decays into tt tt very close to the interac- 
tion point. The Kl decays into further away resulting in four photons which are detected by the 
calorimeter [11]. 

The Ks meson decays into tt+tt" near the collission point. The Kl meson decays into tt'^tt^ and 
then these pions decaying through electromagnetic interaction to 47 quanta. Those four photons 
are detected by the EmC. Thus Kl Stt" — > 67 with two merged clusters may obscure a signal 
from the Kl 27r'' reaction. 
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Implementation of the full barrel geometry into the FLUKA program allowed us to study the 
effect of merging and splitting quantitatively. And hereafter we will present the results for the 
reconstruction of clusters from the rj 37r° and Ks ^ Stt" in view of the merging and splitting 
which estimation is crucial for the study of the 77 7r°77 and Ks 37r° decays. 




active cells distribution 



Figure 5.3: Distribution of the number of active cells for reaction: e e ^ 777 ^ Stt 7 ^ 77. 

Statistics is equal to 10 000 events. 

The Fig. 15.31 shows a distribution of the total number of cells which given signal. Whereas 
distribution of the total number of preclusters and the total number of reconstructed clusters are 
presented in the next Figure. 
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Figure 5.4: Distributions of the reconstructed pre-clusters (left) and total reconstructed clusters (right), 
for reaction: e^e~ — > 777 ^ 37r''7 77. Statistics is equal to 10 000 events. 



The reconstruction possibilities of the clustering algorithm are presented in Fig. 15.41 One can 
see that the most probably the clustering algorithm reconstructs 5 or 6 clusters even though 7 
photons had hit the detector. 



5.2. Splitting of clusters for Ks — ► Stt" as a background for Kg — » Stt" 44 

The frequencies of reconstructed clusters are presented in Table 15.21 
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Table 5.2: Efficiency of the reconstructed pre-clusters and clusters for e+e ^ 777 ^ Stt 7 — > 77 

channel. Statistics equal to 10 000 events. 

One can see that the topology of this channel could be interpreted both with 4 7 quanta 
(15.13%) as well as with 6 photons (23.26%). In the first case this channel could constitute a 
background for r\ — > 7r°77 reaction (merging effect) . 

5.2 Splitting of clusters for Ks 27r° as a background for Ks Stt" 

Splitting is related with the features of the EmC detector and with reconstruction possibilities 
of the KLOE clustering algorithm which not always correctly reconstructs positions and energies 
of particles interacting with detector material and may assign more then one cluster to signals 
induced by a single particle. 

Fig. 15.51 shows that data with 67 quanta in the final state gives in most cases much lower ^ 
when fit under the hypothesis of 27r° then for the 37r'' hypothesis whereas for a real 67 events 
originating from the Ks 37r° decay one expects to be small and xItt to be large as it is seen 
in the sample of the MC Ks -> 37r° data (Fig. [531). 



5.2. Splitting of clusters for Kg 27r" as a background for Kg — > Stt" 



45 



2% 



Data 

MC Ks ^ 371° 



6 ys events !! 




Figure 5.5: Disribution of the variable determined under Stt" and 27r'' hypotheses of the MC Ks Stt" 
signal and of the data [3]. 

Therefore, the Ks 27r'^ 47 decay with two spUtted clusters gives the same event topology as 
expected for the Ks Stt" 67 reaction. We investigated e+e^ cj) ^ KlKs /vL27r" 
Kl^'^ reaction with condition that all four 7 quanta hit the barrel calorimeter. We performed 
approximation that Ks meson decays exactly in the coUision point. However, we suppose that 
this approximation hasn't significant infiuence on the results because a distance between collision 
point and hit position at calorimeter surface is at least 2 meters. We present the distribution of 
the active cells in Fig. 15.61 
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Figure 5.6: Distribution of the number of active cells for e^e ^ — » KlKs — > Kl^tt — > A'_t47 
reaction. Statistics is equal to 10000 events. 

We also investigated the distribution of pre-clusters and clusters for this channel (Fig. 15. 7p . 
One can see that in 18.08% a clustering algorithm reconstructed 67 quanta. Therefore for this 
case the topology of investigated channel with 6 reconstructed photons (two splitted clusters) is 
the same as for Kg ^ 37r° reaction and it constitute a background for this decay. 
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Figure 5.7: Distributions of the number of reconstructed pre-clusters (left panel) and total number of 
reconstructed clusters (right panel), for e+e^ '"''^ r^-o_o 

equal to 10000 events. 



KlKs — > Kl2it — » Kl^j reaction. Statistics is 
The fractions of numbers of reconstructed clusters are presented in the Tab. 15.31 



Number of reconstructed 
preclusters and clusters 


Frequency for 
preclusters (%) 


Frequency for 
clusters (%) 
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0.26 
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10.47 


9.38 
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18.94 


18.08 
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4.84 
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0.94 


0.86 
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0.06 


0.21 


10 


0.01 


0.01 



Table 5.3: Efficiency of the reconstructed pre-clusters and clusters for 
e+e~ ^ — > KlKs Kl2tt" -> A'l47 channel. 

One can see that correct reconstruction of four clusters constitutes only about 26% of total 
number of events. Both merging and splitting effects could be reduced by the better shower shape 
reconstruction induced by a finer read-out granularity. But also the identification of particles based 
on the shower shape analysis would greatly take advantage of an increase in the calorimeter gran- 
uality [3]. 
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5.3 Reconstruction efficiency with the KLOE clustering algorithm 

In order to estimate the reconstruction efficiency of the present KLOE clustering algorithm 
we investigated detector response for two 7 quanta with momentum of 500 MeV/c. The response 
was generated as a function of the distance between the 7 quanta. Subsequently we have re- 
constructed a total number of clusters which were recognized by clustering algorithm. The two 
investigated situations are shown in Fig. 15.81 and in Fig. 15.91 In the first case the distance be- 
tween photons on x axis is cm and the distance on y axis was varried between and 200 cm. 
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Figure 5.8: Hit positions of 7 quanta for AX = cm. Figure 5.9: Hit positions of 7 quanta for 

AY = 200 cm. 

If the reconstruction efficiency of the KLOE algorithm was 100%, this algorithm would always 
reconstruct two clusters which we interpret as two photons which hit the module. But results 
of simulations show (see Fig. 15.111 and Fig. I5.10p that reconstruction efficiency is changing as a 
function of distance between these two particles. 
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Figure 5.10: The efRciency of the two clus- 
ters reconstruction as a function of the distance 
along y axis provided that x coordinate for both 
photons was the same. 



Figure 5.11: The efRciency of the two clusters 
reconstruction as a function of the distance be- 
tween photons along the x axis. 



5.4. Reconstruction efficiency as a function of tlie photocathode quantum efficiency 
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As expected the reconstruction efficiency is increasing with increasing distance between two 
7 quanta. This result shows that the biggest problem with reconstruction is in the case when 
particles hit with the small difference in the x coordinate. 

5.4 Reconstruction efficiency as a function of the photocathode 
quantum efficiency 



We have performed investigations of the influence of the quantum efficiency of the photocatode 
on the merging and splitting effects. We carried out our simulations using a vertex generator for 

37. Then we simulated the calorimeter 



reactions: e+e 



777 87 and e+e 



response using FLUKA program, reconstructed PM response by means of the DIGICLU procedure 
and flnally reconstructed clusters by means of the present KLOE clusterisation algorithm. The 
studies presented in this section were performed for one module only. The efficiency is defined by 
the following formula: 

iV(> Iduster) 



where: N is the tot 
that we take events 



iV(> Ihit) 
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Figure 5.12: Reconstruction efficiency as a function of the energy and quantum efficiency. Courtasy of 
E. Czerwinski [51]. 

The results of simulations for the single calorimeter module are shown in Fig. 15.121 We have 
studied [54] the reconstruction efficiency as a function of the energy in the cluster for five values 
of the quantum efficiency as depicted inside the Figure. QE = 1 represents results for a present 
quantum efficiency of photomultipHers used at KLOE, amounting to 23% [51,60]. Remaining 
values of QE indicates a factors by which the efficiency was increased in our investigations. From 



5.5. Merging and splitting probabilities as a function of the quantum efficiency 
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these studies one can see that as expected with the increasing value of the quantum efficiency for 
photomultipliers the efficiency of reconstruction is increasing too. 

5.5 Merging and splitting probabilities as a function of the quantum 
efficiency 

Subsequently, we determinated the influence of value of the quantum efficiency on the magni- 
tude of the merging and splitting effects. The results of these simulations are shown in Fig. 15.131 
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Figure 5.13: Merging and splitting studies with single EmC module. Rmsmax is a value of the width 
(RMS) of the time distribution of the cells above which the clustering algorithm tries to break a cluster [54]. 
At present KLOE algorithm this value is equal to 2.5 ns. Courtasy of E. Czerwihski [51]. 

The merging probability was defined according to the following formula: 

, , , N(l cluster and > 2 hits) , , 

merg. probability = — ; — — - (5.2) 

^ ^ N(> 1 cluster and > 2 hits) ^ ^ 

and the splitting probability was defined as: 

N(> 2 clusters and 1 hits) 
split, probability = , ,1,,,,, 1 hits) ^''^^ 

Fig. 15.131 shows that the splitting probability increase with increasing of the quantum efficiency of 
photomultipliers [54] but one can see that merging effect is much less sensitive to the changing of 
the quantum efficiency. 
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6. Conclusions 



We have built a full simulation of the Barrel Calorimeter of the KLOE detector with FLUKA 
Monte Carlo package. This implementation of the geometry was used for simulations of the phys- 
ical energy response of the calorimeter. The Monte Carlo cooperates with vertex generator which 
reproduces kinematics of the physical reactions. 

We estimated influence of merging and splitting of clusters for reconstruction possibilities of 
the particles on the whole barrel calorimeter with a present clustering algorithm. In particular 
we estimated these effects as a function of the quantum efficiency of photomultipliers [54]. As a 
results we found that increase of the quantum efficiency does not improve significantly the cluster 
reconstruction in view of merging and splitting effects. 

The prepared program with geometry of the 24 modules of the barrel calorimeter enables fast 
simulations of a choosen reaction channel and facilitates studies of the detector response to differ- 
ent reactions separately. 

This geometry may be used to test response of photomultipliers with higher granularity of cells, 
especially for testing performances of the multianode photomultiplier. 

Implemented geometry of the barrel calorimeter in FLUKA reproduces the merging and the 
splitting effects, and permits also to study the efficiency of clusterisation algorithm at the edges of 
the modules in the barrel calorimeter. These all usefull options come out from the fact that the 
smallest parts of detector (the fibers structure) were realistically implemented in FLUKA Monte 
Carlo. 

We investigated a shapes of peaks for time distribution of reconstruction of clusters, the shapes 
are changing as a function of distance between photons. This information could be very useful for 
the upgrade of the clustering algorithm [61]. 

We estimated also the infiuence of merging and splitting of clusters for e+e' 4> ^ K^Ks 
Kl2'k^ KlAj and e^e^ cj) ^ r]j ^ Stt'^j 7j reactions. The studies were made determin- 
ing the response of the whole barrel calorimeter with DIGICLU program. 

Finally, it is worth mentioning that our study encouraged the authors of the FLUKA to upgrade 
the program and in the near future a new version permitting to replicate the cells more than once 
will be available. 
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A. Kinematic fit procedure 



The kinematic fit procedure is used to correct measured quantities Xi {xi...Xi...Xn) in case of 
redundancy in the measured set of variables. In most cases Xi is normal distributed around xf'"'^ 
with standard deviation ai. 

1 ^ 

^ J^^J2^X, ~ xf (A.l) 

i 

(A.2) 

where a; is a average of the all Xi values. Variables x*'""^ could be constrained by some physical 
law (in our investigations by energy and momentum conservation) in this case they must satisfy k 
equations, which can in general be written in the form: 

F,{xf-^...,x':n^O j = l,-,k (A.3) 

The main purpose of using fit procedure is to find a new approximation (fii) of x*/"^'^ value by 
minimizing the quantity: 

i « 

and imposing the (|A.3p constraints on ^i. This purpose can be achieved using langrangian multi- 
pliers method, that consists of minimizing the quantity: 

= E ^-^^^ + E ^.^.(^1' (A-5) 
' j 

related with the variables fii and Xj [6]. To minimize a value an iterative procedure is needed 
to be used. 

The detailed description of the distribution can be found e.g. in [56]. 
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B. Monte Carlo Methods 



Monte Carlo simulations are methods of solving numerical assignments with the aid of adequate 
constructed statistical computations. They are frequently used in experimental and theoretical 
physic [57]. Although computer memory and processor performance have increased dramatically 
over the last two decades, many physical problems are too complicated to be solved without 
approximations of the physics processes, quite apart from the approximations inherent in any 
numerical method. Therefore, most calculations done in computational physics involve some degree 
of approximation [58] . 

Solving a physical problem often requires to solve an ordinary or partial differential equation. 
This is the case in for example classical mechanics, electrodynamics, quantum mechanics and fluid 
dynamics. On the other hand in statistical physics we must calculate sums or integrals over large 
numbers of degrees of freedom. 

Whatever type of problem we try to solve, it is only in very few cases that analytical solutions 
are possible. In most cases we therefore resort to numerical calculations to obtain useful results [58]. 

These methods are used in physic for the last sixty years [57]. 



55 



56 



C. Energy distributions for 7 quanta 



We have investigateco the distribution of the energy and the position of the 7 quanta on the 
bottom surface of a single calorimeter module. We performed simulations for reactions: e+e^ 
r/7 ^ 37 and e+e^ cj) ^ rj-f ^ 37r°7 77 . The result is presented in Fig. IC.ll in upper 
and lower panel, respectively. 




Figure C.l: Simulated distributions of the hit position at the surface of one calorimeter module (left 
panel) and energy spectra (right panel) of 7 quanta from the reactions e^e^ cj> rij (upper 
panel) and e^e~ —></!> ^ 777 ^ 37r°7 — > 77 (lower panel). 

In the left column in Fig. IC.ll one can see hit distributions for particles entering the module. 
In the right panels in this figure the sharp signals orginate from a monenergetic 7 quanta from (j> 
meson decay with energy equal to 363 MeV. In the case of the e+e~ ^ 777 ^ 37r"7 77 

^We used the BXDRAW subroutine in this investigations [32]. 
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reaction we observe more events at the lower energy range. These signals orginate from 7 quanta 
from 7r° meson decays. As a further example in Fig. IC.2I we present energy spectra separately for 
the radiative photon and for the photon from the 77 and pion decays. 
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Figure C.2: The energy distribution for photons from reactions: e+e ^ <^ ^ 7^7 ^ 87 and 

e^e~ ^ <j) ^ rj^ ^ 37r"7 — > 77. 

In the left panel the monoenergetic signal for photon from (j) meson decay is shown. The middle 
panel presents energy distribution for photons from the 77 meson decay and the third panel presents 
the energy distribution for 7 quanta from 7r° mesons decays. 



D. Time distribution for single and mult i- gamma hits 



In this appendix we investigated distributions of the sum and difference of times of Hght signals 
from the sides of the module. Our studies show that distributions are changing as a function of 
the distance between 7 quanta hitting the module. This information could be used to improve 
efficiency of the KLOE clustering algorithm in the near future [61]. In Fig. ID. II result for the 
situation where one 7 quanta hit in the middle of the calorimeter module is presented. One can 
see only one peak for both difference and sum of time of signals from sides A and B. 
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Figure D.l: Distribution of sum and difference of time for one 7 quanta. Ta and Tg denote the time 
where a signal from scintillating fibers comes into a photomultipliers on side A and side B, respectively. 

For the two 7 quanta which hit the surfacil] at the same x but different y the result is different 
(see Fig. ID.2p . One can see that distribution of difference of times for signals from side A and 
B consists of three separated peaks and the distribution of the sum of times from sides A and B 
shows two maxima. 



^Two studing cases of hit positions of photons entering the calorimeter surface are presented in Fig. 15.81 and 
Fig. 15.91 Positions of 7 quanta are marked as black circles. 
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Figure D.2: Distribution of sum and difference of time for two 7 at the same x position but different y. 

Next we investigated the situation where distance on x axis between two 7 quanta is equal to 
15 cm. In this case we observe only one maximum in the distribution of the sum and two maxima 
in the distribution of the difference of times (Fig. ID.31) . 
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Figure D.3: Distribution of sum and difference of time for two 7 hitting the module with a distance of 
15 cm along the x axis. 



In the following we present an example of results of cell distributions simulated for several 
distances on x axis between places where photons hit the module and for several angle between a 
surface and a photon direction. The results of the sum of energy deposited in cells for side A and 
B we present using a displayer progranjf]. Fig. ID. 4l presents a photomultipliers energetic response 
for one 7 quanta which hit the middle of the module at an angle of 90 degrees. 



^This program written in a Fortran code [54] is based on the CERN library tool [62]. 
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Figure D.5: Distribution of cells for two 7 

Figure D.4: Distribution of cells for one 7 

quanta hitting with 90 degrees with a distance 

quantum hitting with 90 degrees in the middle 

on y axis equal to 200 cm and with the same 

of the module. 

coordinate on x axis. 



Comparing this result with the case (see Fig. ID.SP where two 7 quanta hit the module at a 
distance on y axis equal to 200 cm but with the same x coordinate, we can see that it is not 
possible to distinguish these two situations due to the fact that the second case can be interpreted 
for example as one particle which hit the module but with higher energy. However, we can easily 
identify two particle events when a distance on the x axis between hit positions is larger than 15 cm 
(see also Fig. 15.11]) . In Fig. ID.7l a distribution for the case when the distance between two photons 
was equal to 25 cm is presented. 




2 4 6 8 10 12 2 4 6 8 10 12 



COLUMNS COLUMNS 



Figure D.6: The distribution of cells for one Figure D.7: The distribution of cells for two 
7 quantum hitting the module with the angle 7 quanta hitting the module with 25 cm dis- 
equal to 30 degrees with respect to the surface. tance along x axis. 

We simulated also a situation where photon hit the module with the angle equal to 30 degrees 
(see Fig. ID.6[) with respect to the surface. One can see two cells on the right side on this plot 
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which aren't connected to the main cluster. It is worth to mention that such distributions may be 
interpreted as two hits. 



E. Estimation of probability for mult i- gamma hits at a single 
calorimeter module 



Using FLUKA Monte Carlo with vertex generator we were able to estimate a probability on 
the multi- gamma hits at calorimeter module orginating from the rj meson decays. We estimated 
a multiplicity of the 7 quanta on one module of the KLOE calorimeter. By multiplicity we define 
the number of 7 quanta hitting one module (Fig. lE.l]) . 
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Figure E.l: The multiplicity of the 7 quanta on a single calorimeter module. For reactions described in 
the figure. 



If it is equal to zero, it signifies that our tested module did not seen any 7, if it is equal to one 
than the module registered one gamma quantum. Due to the axial symmetry of KLOE this result 
is valid for all detection modules. Thus, we see that for the ^ 777 87 reaction only for 0.6% 
of the total number of events the 27 quanta hit one module and in the rest 99.4% only one 7 is 
registered by one module. However, for the reaction 4> 777 7j where 77 meson decays into 37r° 
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the situation changes drastically and already in 6.3% cases there is two or more 7 quanta hitting 
the module. 

In the situation where more than one photon hit the calorimeter module it is possible that 
merging effect will appear for some distances between hit positions of particles on the surface. The 
estimations of the order of multiplicity effect was shown in Fig. lS.lll and Fig. 15.101 In the following 
we will describe the distribution of the distance between the 7 quanta in the case when two 7 
hits the module. This result is relevant in view of the reconstruction possibilities of the clustering 
algorithm (see Fig. 15.111 and Fig. I5.1Q|) . The result for the rj 37r° is shown in Fig. IE. 21 The 
corresponding reactions are written inside the figures. 
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Figure E.2: Distance between two gamma quanta from 77 — > Stt in one calorimeter module. 

And the result for the 77' 37r° is shown in Fig. IE. 31 The distribution determined for the 77' is 
not very intuitive. The maximum at 1.7 m is due to the fact that the pions from the 77/ decay are 
very fast and both 7 quanta from one pion can hit one module. The pion decays immediately in 
the target and the distance between the two 7 quanta from its decay at the calorimeter surface is 
indeed around 1.7 m. 
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Figure E.3: Distance between two 7 quanta from reaction 77' Stt" on one calorimeter module. 
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Figure E.4: Distance between two 7 quanta on one calorimeter module for reactions as depicted in the 
figures assuming that the "pion" mass is equal to 10 MeV. 

We did also one more check by decreasing in the simulations the mass of the pion to 10 MeV. 
In such a case these "pions" are much faster and the laboratory angle between gamma quanta 
is decreasing significantly, and the distance between them in the calorimeter is as expected to be 
much smaller. The result of this simulation is shown in Fig. IE. 41 
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F. Example of the event reconstruction for 

e+e^ ^ ^ ?77 ^ 37r^7 77 reaction 

The capabilites of the prepared program to reproduce a merging and spHtting effects in light 
output from scintilating fibers are ilustrated in Fig. IF.ll 
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Figure F.l: Splitting of energy deposits in the barrel calorimeter. 

One can see seven significant energy deposits belonging to 7 photons from the e+e^ (f) ^ 
rjj 37r"7 77 reaction, and also some small deposits of energy which don't orginate directly 
from photons are well seen. The splitted energy deposits may cause reconstruction of the false clus- 
ters. On the other hand two very close tracks of photons (near z = -250 cm) may be reconstructed 
as one merged cluster since the distance between particles is only about 4 cm, leading with a large 
probability to false reconstruction (see Fig. IF.1|) . After simulating the energy response of the scin- 
tillating fibers, on 24 modules with FLUKA, we extended appropriately the DIGICLU program 
in order to perform reconstruction of signals from photomultipliers on the whole barrel calorimeter. 
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Figure F.2: Reconstruction of clusters for reaction: e^e — > <^ ^ 777 — > 37r'^7 77. Statistics is equal to 
1000 events. 

Fig. IF. 21 presents distribution of reconstructed clusters on the barrel part of the calorimeter 
for the 1000 e+e~ ^ 777 — > 37r°7 — > 77 events simulated with vertex generator and FLUKA 
program. One can see that the clusters aren't reconstructed on the edges of modules. 



G. Energy deposition as a function of azimuthal angle 



We studied energy deposits in scintillating fibers as a function of the azimuthal angle of 7 
quanta producted in the interaction region. The results for the upper part of barrel calorimeter is 
shown in Fig. IG.ll 
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Figure G.l: Energy deposits on the upper half of the barrel calorimeter as a function of cylindrical angle. 




Figure G.2: Schematic view of the 4> angle distribution. 



One can see that at the edges of the module a maximum of the energy deposits appears. To 
some extent it may be explained by the fact that on the edges the value of the solid angle for the 
same 4> angle interval is larger than in the middle (see Fig. IG.2( ). 
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Figure G.3: Energy deposits on the barrel calorimeter as function of the cylindrical angle in the range 
from to 22.5 degree. 

In the simulations we did not include fibers which were near the edge since in a real detector the 
trapezoid shape of the module was acompHshed by cutting some parts from rectangular piece [11] 
and hence the scintillating fibers at the sides were broken. In Fig.lOSl one can see that for value of 
an angle around 0.15 rad particles didn't deposited energy. This is because this position is exactly 
between the edge of the first and second trapezoid modules. 



H. Definition of the coordinate system 

The coordinate systems which were used in simulations are shown in Fig. IE. II The blue axis 
denote the frame which was used for simulations with FLUKA and with the VERTEX GENERA- 
TOR. The red coordinate system which is the same as coordinate frame of the KLOE detector has 
been used in simulations of clusterization effects and simulations of the photomultipliers response 
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Figure H.l: The coordinate systems used in simulations. 



Throughout this thesis we plot all figures using the coordinate frame of the FLUKA. 



^Digiclu frame is called the reconstruction geometry system. 
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The world is given to me only once, not one ex- 
isting and one perceived. Subject and object are 
only one. The barrier between them cannot be 
said to have broken down as a result of recent ex- 
perience in the physical sciences, for this barrier 
does not exist. 



Swiat dany mi jest tylko raz, nie ten istniej§,cy 
i nie ten doswiadczany. Podmiot i obiekt s§, 
jednym. Bariera pomigdzy nimi nie moze bye 
zlamana jako wynik ostatnich doswiadczeh w 
naukach fizycznych, bariera ta nie istnieje. 



Erwin Schrodinger (1887 - 1961) 



Education is what remains after one has forgotten 
everything he learned in school. 



Edukacja jest tym co zostaje jak zapomni si§ 
wszystko czego nauczylo si§ w szkole. 



Albert Einstein (1879 - 1955) 
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